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Abstract 

The present work extends our short communication [lj. For smooth 
marginally outer trapped surfaces (MOTS) in a smooth spacetime we de- 
fine stability with respect to variations along arbitrary vectors v normal 
to the MOTS. After giving some introductory material about linear non 
self-adjoint elliptic operators, we introduce the stability operator L v and 
we characterize stable MOTS in terms of sign conditions on the princi- 
pal eigenvalue of L v . The main result shows that given a strictly stable 
MOTS So C So in a spacetime with a reference foliation S^, there is an 
open marginally outer trapped tube (MOTT), adapted to the reference 
foliation, which contains Sq. We give conditions under which the MOTT 
can be completed. Finally, we show that under standard energy conditions 
on the spacetime, the MOTT must be either locally achronal, spacelike or 
null. 



1 Introduction 

The singularity theorems of Hawking and Penrose assert the presence of incom- 
plete geodesies in the time evolution of Cauchy data with physically reasonable 
matter containing a trapped surface [2]. Studying the structure of these singu- 
larities and understanding whether generically they entail a loss of predictability 
power of the theory have become central issues in Classical General Relativity. 
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More specifically, the weak cosmic censorship conjecture asserts, roughly speak- 
ing, that in the asymptotically flat case the singularity will be generically hidden 
from infinity by an event horizon and that a black hole will form. Since this 
conjecture aims precisely at showing that a black hole forms, any sensible ap- 
proach to its proof should not make strong a priori assumptions on the global 
structure of the spacetime. It is therefore necessary to replace to concept of black 
hole, which requires full knowledge of the future evolution of a spacetime, with 
a quasi-local concept that captures its main features and that can be used as a 
tool to show the existence of a black hole. 

This approach has been successfully applied in spherically symmetric space- 
times with matter fields, where the existence of a complete future null infinity 
and an event horizon can in fact be inferred from the presence of at least one 
trapped surface in the data [3], plus some extra assumptions. An important tool 
in the analysis is to study the sequence of the marginally trapped surfaces bound- 
ing the region with trapped surfaces within each slice of a spherically symmetric 
foliation. 

Quasi-local versions of black holes are important not only in the context of 
cosmic censorship, they are relevant in any physical situation involving black 
holes where no global knowledge of the spacetime is available. An outstanding 
example is the dynamics and evolution of black holes. In the strong field regime, 
these evolutions are so complex that they can only be approached with the aid 
of numerical methods. In most cases, numerical computations can only evolve 
the spacetime a finite amount of time, which makes the global definition of black 
hole of little practical use. A quasilocal definition becomes necessary even to 
define what is understood by a black hole in this context. More importantly, 
such a definition is crucial in order to be able to track the location of the black 
holes and to extract relevant physical information from their evolution. Over 
the last few years, marginally outer trapped surfaces and the hypersurfaces in 
spacetime which they sweep out during a time evolution, have become standard 
as quasilocal replacements of black holes and have been studied extensively, using 
both numerical methods (see, e.g. [U El E]) as well as analytically, with either 
mathematical [7J or more physical scope [8] . See [9] for a review of some of these 
issues. However, many open problems still remain. Only in spherical symmetry a 
rather complete picture has been obtained thanks to [3j [TOj [TT] . In general even 
finding examples is not easy, see however [TJl EH EE] for interesting non-spherical 
examples. 

A marginally outer trapped surface (MOTS) is a spacelike surface of codi- 
mension two, such that the null expansion 6 with respect to the outgoing null 
normal vanishes. The notions of "outer" and "outgoing" are simply defined by 
the choice of a null section in the two dimensional normal bundle of the surface. 
We call a 3-surface foliated by MOTS a marginally outer trapped tube (MOTT) 
[TSl [TB]; an alternative terminology adopted in [Ij is "trapping horizon", c.f. [17]. 

The wide range of applicability of MOTS motivated us to study their propa- 
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gation in spacetime from an analytical point of view and in a general context, i.e. 
assuming neither symmetries nor the presence of any trapped regions a priori. In 
the context of the initial value problem in general relativity, namely in a smooth 
spacetime foliated by smooth hypersurfaces S t , it is natural to ask the following: 
Given a MOTS So on some initial leaf E , does it "propagate" to the adjacent 
leaves S 4 of the foliation? In other words, is there a marginally outer trapped 
tube starting at So whose marginally outer trapped leaves lie in the time slices 

It turns out that the key property of a MOTS So relevant for this question 
is its "stability" with respect to the initial leaf So- This concept has interesting 
applications even when considered purely inside a hypersurface S , in particular 
for the topology of So, and also for the property of being a "barrier" for weakly 
outer trapped and weakly outer untrapped surfaces (defined by 9 < and 9 > 0, 
respectively). We shall discuss these two issues, which both originate in the 
work of Hawking [T51 US], before turning to the question of propagation of S 
off S . Hawking's analysis was extended by Newman [20] who calculated the 
general variation 5 V of the expansion 9 with respect to any transversal direction 
v. A central issue in these papers was to show that stable MOTS have spherical 
topology in the generic case. The classification of the "rigidity case" , in which the 
torus is allowed, was investigated first for minimal surfaces [2T] and subsequently 
also for generic MOTS and in higher dimensions [22} [231 El] ; see [25] for a review. 

A key tool, both for the topological issues as well as for the present purposes, 
is the linear elliptic stability operator L v defined by 5^ v 9 = L v ip for MOTS (in- 
troduced in [TJ, in deformation form already present in [2QI [22]) where v is a 
suitably scaled vector. L v is not self-adjoint in general, except in special cases as 
for example when the MOTS lies in a time symmetric slice. Nevertheless, linear 
elliptic operators always have a real principal eigenvalue, and the corresponding 
principal eigenfunction can be chosen to be positive [26], [271 12H [22] . While we 
define stability and strict stability of MOTS in terms of sign conditions for pre- 
ferred variations 5 V 9, we can show that this definition is equivalent to requiring 
that the principal eigenvalue of L v is non-negative or positive, respectively. In 
[I] we also showed that strictly stable MOTS are barriers for all weakly outer 
trapped and weakly outer untrapped surfaces in a neighbourhood. In addition 
to the properties of the stability operator, this result uses a representation of 
MOTS as graphs over some reference 2-surfaces. In terms of this representation, 
the condition of vanishing expansion 9 characterizing a MOTS becomes a quasi- 
linear elliptic equation for the graph function, for which a maximum principle 
holds. A similar application of the maximum principle for the functional 9 is 
contained in the uniqueness results of Ashtekar and Galloway [7J where null hy- 
persurfaces through a given MOTS and their intersection with a given spacelike 
MOTT were considered. 

A barrier property of trapped and marginally trapped surfaces, which com- 
plements the one discussed above has been considered by Kriele and Hayward 
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[30] . They showed that the boundary of the trapped region, i.e. the set of all 
points in a spacelike hypersurface contained in a bounding, trapped surface is a 
MOTS, under the assumption that it is piecewise smooth. By p], this MOTS is 
necessarily stable. Andersson and Metzger [31] recently showed that the bound- 
ary of the trapped region is a smooth, embedded, stable MOTS, without any 
additional smoothness assumption. 

Turning now to our main problem, namely the propagation of the MOTS 
So off S to adjacent slices, we can prove this for some open time interval if 
the MOTS are strictly stable (c.f. Theorem 19. 2D . As a tool we first extend the 
graph representation of So C So to 2-surfaces S t C E 4 . The linearization of 
the expansion operator 9 is precisely the stability operator, and strict stability 
guarantees that this operator is invertible. As So is a solution of the equation 
9 = on S , we could apply in our earlier paper [1] the implicit function theorem 
to get solutions of 9 = in a neighbourhood. Here we cut this procedure short 
by using standard results on perturbations of differential operators [321 E3] whose 
linearizations are elliptic and invertible. Naturally, these results also make use 
of the implicit function theorem. As an easy corollary to Theorem I9.2[ we find 
that the MOTT constructed in this interval is nowhere tangent to the S t . Much 
more subtle results are Theorem 19.41 and Corollary 19.51 Under some genericity 
condition they show in particular that, if St converge smoothly to a limiting 
MOTS S T C S T whose principal eigenvalue A T vanishes, the resulting MOTT is 
everywhere tangent to E T . 

We wish to stress that the existence result contained in [Tj is local in time. The 
attempt to formulate a global result in [Tj required the assumptions (implicit in 
the definitions of that paper) that the MOTS remain compact, embedded, smooth 
and strictly stable during the evolution. In the present work we allow immersed 
rather than embedded MOTS in our existence theorems for MOTT, but we have 
to deal with the other potential pathologies. To do so we apply a recent result by 
Andersson and Metzger [31] which shows that, in four dimensional spacetimes, a 
sequence of stable and smooth MOTS which lie in a compact set such that area 
of the sequence stays bounded, converges to a smooth and stable MOTS. This 
leads to Proposition 19.61 Moreover, if the dominant energy condition holds, it 
is easy to show that the area stays in fact bounded provided the MOTS remain 
strictly stable in the limit. This gives Theorem 19.81 as a sharper version of our 
existence result. 

This paper is organized as follows. In section [2] we explain the most important 
items of our notation. In section [3] we discuss the variation of the expansion, and 
introduce the stability operator. The somewhat technical computation of the 
variation, which simplifies the derivation by Newman, is given in appendix |A] 

We proceed with some technical material on linear elliptic operators with 
first order term, cf. section |H Here and in appendix [B] we give an exposition 
of, in particular, the Krein-Rutman theorem [35] on the principal eigenvalue and 
eigenfunction of linear elliptic operators, and the maximum principle for operators 
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with non-negative principal eigenvalue |36j. We continue in section[5]by discussing 
in detail stability definitions for MOTS, in particular the relation between the 
variational definitions and the sign condition on the principal eigenvalue, and we 
give a result on the dependence of stability on the direction. Sect. |6] contains the 
graph representation of a MOTS. In Sect. [7] we describe the barrier properties of 
MOTS which satisfy suitable stability conditions, along the lines sketched above, 
slightly extending our earlier paper p] . In section [HJ we show, roughly speaking, 
that strictly stable MOTS inherit the symmetry of the ambient geometry, and 
that the same is true for the principal eigenfunction. In the final section [9] we 
prove existence for MOTT in the three theorems 19.21 I9T41 and 19.81 already sketched 
above. Our final Theorem 19.91 is a slight extension of a result of pQ and shows 
that under standard energy conditions, suitable (non-) degeneracy conditions for 
the initial MOTS Sq and for spacelike or null reference foliations, the MOTT 
through Sq is either spacelike or null everywhere on Sq. 

2 Some basic definitions 

A spacetime (Ai,g) is an n-dimensional oriented and time-oriented Hausdorff 
manifold endowed with a smooth metric of Lorentzian signature +2. Some results 
below require n = 4. S will denote an orientable, closed (i.e. compact without 
boundary) codimension 2, immersed submanifold of M. with positive definite first 
fundamental form h. An object with all these properties is simply called "surface" 
throughout this paper. The area of S will be denoted by Spacetime tensors 
will carry Greek indices and tensors in S will carry capital Latin indices. Our 
conventions for the second fundamental form (-vector) and the mean curvature 
(-vector) are K(X,Y) = -(V^y) 1 and H = tihK. Here X, Y are tangent 
vectors to S, _L denotes the component normal to S and V is the Levi-Civita 
connection of g. S will always be assumed to be smooth unless otherwise stated. 

The normal bundle of S is a Lorentzian vector bundle which admits a null 
basis {l a , k a } which we always take future directed, smooth and normalized so 
that l a k a = —2. This basis is defined up to a boost l a — > nl a , k a — > k _1 /c q , 
k > 0. The null expansions of S are Q\ = H a l a and 6^ = H a k a and the mean 
curvature in this null basis reads H a = — | (6kl a + 6ik a ). 

Definition 2.1. A surface S is a marginally outer trapped surface (MOTS) 
if H a is proportional to one of the elements of the null basis of its normal bundle. 

This condition is more restrictive than just demanding H a to be null because 
it excludes the possibility that H a points along l a in some open set and along 
k a on its complement (c.f. [37]). The null vector to which H a is proportional 
is called l a , and the direction to which it points is called the outer direction (in 
the case H a = 0, both l a and k a are outer directions). In other words, the term 
outer does not refer to a direction singled out a priori, but to the fact that we 
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only have information about, or we are only interested in, one of the expansions. 
Equivalently, S is a MOTS iff $i = 0. If furthermore H a is either future or past 
directed everywhere S is called marginally trapped. Hence a marginally trapped 
surface satisfies Q\ = and either 9k < or 9k > everywhere. Next, S is called 
weakly outer trapped iff at least one of the expansions in non-positive, say 9i < 0. 
Weakly outer untrapped surfaces satisfy the reverse inequality. Finally, in order 
for S to be a future (past) trapped surface we require that the strict inequalities 
Q\ < and 9k < 0, (9i > 0, 9k > 0) hold. Since we will only deal with 9% from 
now on, we simplify the notation and refer to it simply as 9. 

A marginally outer trapped tube Q is a smooth collection of MOTS. More 
precisely, we state the following 

Definition 2.2. Let I C R be an interval. A hypersurface Q, possibly with 
boundary, is a marginally outer trapped tube (MOTT) if there is a smooth 
immersion $ : S x I — > M , such that Q — $(<S X I) and 

(i) for fixed s <E I , $(5, s) is a MOTS with respect to a smooth field of null 
normals l a on Q and 

(ii) $*(<9 S ) is nowhere zero. 

Suppose (M,g) contains a foliation by hypersurfaces {S t } te j. A MOTT Q is 
said to be adapted to the reference foliation {St} if for each s G / , it holds 
that S a [ s ) = $(<S, s) is a MOTS in E (T ( a ) ; where a : / — » J is a smooth, strictly 
monotone function. 

Remark. Note that if / contains at least one of its boundary points, then the 
MOTT Q is a hypersurface with boundary. 



Condition (ii) is required in order to allow self-intersections of the MOTS but 
not in the direction of propagation. For embedded MOTS, the MOTT is also 
embedded and its definition is equivalent to Hayward's "trapping horizons" [18J. 
The terms "dynamical horizon" and "isolated horizon" [15l [16] are particular 
cases in which the causal structure is restricted a priori . 



3 Varying the expansion 

A fundamental ingredient in our existence theorem is the first order variation 
of the vanishing null expansion 9 of a MOTS. This variation was given in full 
generality by Newman in [20J for arbitrary immersed surfaces (not necessarily 
MOTS). We give here a simplified derivation. We first have to introduce some 
notation. 

Let V„ and G a p denote the covariant derivative and Einstein tensor of (A4,g) 
respectively. Let («S, h) be a spacelike codimension-two surface immersed in M. 
(not necessarily closed), e a A any basis of the tangent space of S, Da the covariant 
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derivative on (S, h) and Rs its curvature scalar. Let us fix a null basis {l a , k a } 
in the normal bundle of S and define a one-form on S as 

SA = ~^k a W e A l . 

We shall calculate how 9 changes when the surface S is varied arbitrarily. This 
variation is defined by an arbitrary spacetime vector q a defined along S. More 
specifically, let G I C K. be an open interval and $:5x/^jVlbea 
differentiable map such that for fixed a, $(•, a) is an immersion and for fixed p, 
Xp(a) = $ a (p, a) is a curve starting at p 6 5 with tangent vector q a (p). Define 
the family of two-surfaces S a = $(<S, a). Let be a nowhere zero null vector on 
the normal bundle of S a which is differentiable with respect to a, and let 9 a be the 
null expansion of each surface S a . The variation of 9 a , defined as S q 9 = d cr 9 cr \ r7= o, 
depends only on q a and on the null vector field l a and its first variation (if S is 
a MOTS this last dependence also drops out), but not on the details of the map 
$. For a MOTS, the variation is moreover linear in the sense that 

$aqi+q 2 9 = aS qi^ + &q-fi (!) 

for any constant a, while in general 5^9 ^ ipo~ q 9 for functions ip. 

It should be noted that in the context of trapping and dynamical horizons, 
derivatives of 9 have been employed frequently (for instance in the definition of 
outer trapping horizon by Hayward p2] or in the uniqueness results by Ashtekar 
and Galloway [7]). These are not the variations we are considering in this paper. 
The former are derivatives of a scalar function defined in a neighbourhood of the 
horizon by extending 9 off the horizon, using the Raychaudhuri equation. In this 
case, the derivative is obviously linear with respect to multiplication by scalar 
functions, unlike the geometric variation employed here. 

The variation vector q a can be decomposed into a tangential part g" a and an 
orthogonal part q^ 06 with respect to S, i.e. q a = q ±a +q" a . The normal component 
can in turn be decomposed in terms of the null basis as q ±a = bl a — ^k a where 
b, u are functions on S. The following result, whose proof we have shifted to 
appendix [A], gives the explicit expression for the variation of 9 along q a . 

Lemma 3.1. Let S, l a , 9 and q a = + bl a - \k a be as before. Then, the 
variation of 9 along q is 

S q 9 = a9 + qW{9)-b (K» B K uAB l^l u + G^Fl v ) - A s u + 2s A D A u 
+ ^(R S -H 2 ~ G^k v - 2s A s A + 2D A s A ) , 

where A$ = D A D A is the Laplacian on (S, h) and a = ~^k a d a l^\ a= Q. 

The decomposition of q ±a in the null basis {l a , k a } is natural for the surface S 
as a co dimension-two submanifold in spacetime, and Lemma [3TT1 gives the general 
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variation with respect to arbitrary vectors on S. However, we will later refer the 
variations of S to some foliation of A4 by hypersurfaces E t with S C So ( c -f- 
section Ej), and for this we will employ a natural alternative decomposition of q ±a 
adapted to the foliation. We consider an arbitrary normal vector field v a to S 
which is, at each point, linearly independent of l a , and we impose the normaliza- 
tion v a l a = 1 which does not restrict the causal character of v a anywhere on S. 
v a is uniquely defined by a scalar function V on S according to 

v a = - X -k a + VT. (2) 

We use {v a , l a } as a basis in the normal space. Inverting ([2]) we get k a = 
2{Vl a — v a ). We next define a vector u a = ^k a + Vl a , which is orthogonal to v a 
and satisfies u a u a = —v a v a = —2V. The quantities 

W = K» B K» A %l v + G^m u , (3) 
Y = VK^K^l^ + G^u v (4) 

will appear frequently below. Clearly W is non-negative provided the null energy 
condition holds and Y is non-negative if u a is causal and the dominant energy 
condition holds. 

The following definition introduces an object which plays a key role in this 
paper. 

Definition 3.1. The stability operator is defined by 

L v ij = -A s 4> + 2s A D A ij + Qifc -Y- s A s A + D A s A ^j i/>. (5) 

The following lemma is a trivial specialization of Lemma 13.11 

Lemma 3.2. Let S be a MOTS, i.e. 9 = 0. The variation of the expansion 
9 on S with respect to the null vector ipl a , and with respect to any vector ipv a 
orthogonal to S with l a v a = 1, respectively, are given by 

5^9 = (6) 
5rp v 6 = L v if), (7) 

We note that ([6]) is the Raychaudhuri equation. For arbitrary vectors w a 
orthogonal to S and linearly independent of l a , not necessarily normalized to 
satisfy w a l a = 1, we can define another elliptic operator L w ip = 5^ w 9. In terms 
of the normalized vector v = F~ x w, where F = w a l a we obviously have 

L w i) = L Fv ^ = L v (Fip). (8) 
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Hence, L v and Lp v contain essentially the same information. To see the depen- 
dence of L v on the vector v, i.e. on the function V we calculate, from linearity 
O, and ©, 

L v = l -L_ k - VW. (9) 

Due to the presence of the first order term in ([5]), L v is not self-adjoint in 
general. However, self-adjoint extensions still exist in special cases. For example, 
if sa is a gradient, i.e. sa = Da( for some (, L v is self-adjoint with respect to a 
suitable measure depending on (, c.f. section HJ 

Since L v is linear and elliptic, it has discrete eigenvalues and the corresponding 
eigenfunctions are regular. However, in general, the eigenvalues and eigenfunc- 
tions are complex. Nevertheless, the principal eigenvalue, i.e. the eigenvalue 
with smallest real part, and its corresponding eigenfunction behave in the same 
manner as for self-adjoint operators. In particular, they can be used to give a 
very efficient reformulation of the maximum principle. In the next section we 
collect some material on linear elliptic operators which will be the key tools in 
the subsequent discussion of stability. 

4 Properties of linear elliptic operators 

The results of this section hold for connected compact smooth manifolds S with- 
out boundary, and for arbitrary smooth, linear, second order, elliptic operators 
on S, which can be written as 

L = -A h + 2t A D A + c (10) 

where A h ip = Da (h AB D B ip) , h AB is positive definite and smooth, Da is the 
corresponding Levi-Civita covariant derivative and t A and c are smooth. 

Lemma 4.1. Let L be an elliptic operator of the form ( TiOj) on a compact mani- 
fold. Then, the following holds. 

(i) There is a real eigenvalue \, called the principal eigenvalue, such that for 
any other eigenvalue fi the inequality Re(fi) > A holds. The corresponding 
eigenfunction <fi, Lip = A0 is unique up to a multiplicative constant and can 
be chosen to be real and everywhere positive. 

(ii) The adjoint L^ (with respect to the L 2 inner product) has the same principal 
eigenvalue A as L. 

Applications of Lemma I4~T1 to the stability operator ([5]) in a spacetime will be 
described in the next section. Note that our terminology follows Evans [29]; in 
particular, regarding the sign of L it is opposite to other references [26l [36] cited 
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below. This entails that when comparing with these references one always has to 
interchange "sup" and "inf" when acting on expressions containing L. 

The existence of the principal eigenvalue and eigenfunction was first proven 
by Donsker and Varadhan [26] using parabolic theory and the Krein-Rutman 
theorem [35]. For uniqueness of the principal eigenfunction, c.f. Berestycki, 
Nirenberg and Varadhan [36]. All these papers actually deal with the Dirichlet 
problem for bounded domains in W 1 . However, the proof is easily adapted to the 
case of compact connected manifolds without boundary (it is in fact simpler). For 
completeness and since this result is not widely known, we provide in appendix iBl 
a sketch of the proof in the closed manifold case. The sketch follows the argument 
in Smoller [28] and Evans [29] . 

For self-adjoint operators L = —A h + c, the principal eigenvalue A is given 
by the Rayleigh-Ritz formula 

Ao = inf (u, L u) = inf / (D A uD A u + cu 2 ) r)s, (11) 
U u J s 

where ?]$ is the surface element on (S, h) and the infimum is taken over smooth 
functions u on S with ||m||l2 = 1. 

For non- self-adjoint operators such a characterization is no longer true. How- 
ever, Donsker and Varadhan [26] have given alternative variational representa- 
tions of the eigenvalue, namely 

A = inf sup / if>~ l Lif)^i S) (12) 
ip Js 

A = sup inf ^~ 1 (x)L?/'(x). (13) 

Here the infimum in ([121) is taken over all probability measures fig on S, while 
the suprema are over all smooth, positive functions if) on S. 

To get (H3J from $T2& we first note that the "inf" and "sup" in flT2j) can 
be interchanged (which is non-trivial but follows from a min-max theorem of 
Sion [38], c.f. [26]). This done, the infimum of the integral with respect to all 
probability measures is achieved by a Dirac delta measure concentrated at the 
location where the integrand takes its infimum. 

In order to approach a characterization closer to a Rayleigh-Ritz expression, 
we note that, since probability measures can be approximated by smooth positive 
functions, we can assume /is = u 2 rjs with smooth u > and | \u\ \^ = 1- Following 
Donsker and Varadhan, if) can be decomposed as if) = ue^ '. Direct substitution 
in (TT2"j) and rearrangement gives 

A = inf sup / (D A uD A u + (c + t A t A - D A t A ) u 2 - (D A u + t A ) 2 u 2 ) r] S . 
u w Js 

To reformulate this expression we use the Hodge decomposition t A = D A ( + z A , 
where ( is a function and z A is divergence-free. This decomposition is unique 
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except for a constant additive term in (. The supremum over u only affects the 
last term, and it only depends on z A and not on (. Thus, we need to deter- 
mine inf^ j s (D A u + za) 2 u 2 i]s, for each u. A standard argument shows that the 
infimum is achieved and is given by the solution of the linear elliptic equation 

- D A D A u[u] - 2u~ 2 D a uj[u}D a u = 2 U - 2 z A D A u, (14) 

where we write u[u] to emphasize that the solution depends on u. A trivial 
Fredholm argument shows that this equation has a unique solution satisfying 

J u 2 u[u]rj s = 0. (15) 

It therefore follows that the Donsker-Varadhan characterization of the principal 
eigenvalue can be rewritten as 

A = inf / (D A uD A u + Qu 2 - {D A u[u\ + z A f u 2 ) rj s , (16) 
u Js 

where Q = c + t A t A — D A t and the infimum is taken over smooth functions of 
unit L 2 norm. The symmetrized operator L s = —Ah + Q is self-adjoint and has 
a principal eigenvalue A s given by the Rayleigh-Ritz formula, as in ( fTTI) . Since 
the last term in ( fl6|) is non-positive, the inequality 

A s > A (17) 

follows immediately. This inequality has recently been demonstrated by Galloway 
and Schoen [23] using direct estimates. It is interesting that (fT7|) is a trivial 
consequence of the Rayleigh-Ritz type formula ffl6|) for the principal eigenvalue. 

As a second application of the characterization ( fl~6j) . we note that the last 
term can be rewritten as 

- (D A u[u] + z A f u 2 = (D a lo[u]D a uj[v] - z A z A ) u 2 -2D A (cu[u]u 2 (z A + D a uj[u})) 

after using the equation (Till) and the fact that z A is divergence-free. Integration 
on S gives the alternative representation for A 




Since the last term is now non-negative, dropping it decreases the integrand and 
therefore also the infimum. Thus, defining the alternative symmetrized operator 
L z = L s — zaz a and denoting by \ z the corresponding principal eigenvalue, it 
follows that 

A > A 2 . 

Therefore the principal eigenvalue of a non self-adjoint operator is always sand- 
wiched between the principal eigenvalues of two canonical symmetrized elliptic 
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operators, and we also note that mfs(zAZ A ) < X s — \ z < swp s (zAZ A ) . Obviously, 
when L is self-adjoint (w.r.t. to the L 2 norm) i.e. = 0, the two symmetrized 
operators L s and L z coincide with it. More generally, when t& is a gradient (so 
that L is self-adjoint with a suitable measure) the characterization of A given by 
(TIB]) coincides with the Rayleigh-Ritz expression because whenever z A = 0, the 
unique solution to ffl4|) and ffT5l) is just u[u) = for all u. 

An important tool in the analysis of the properties of the stability operator 
will be the maximum principle. The textbook formulations normally require that 
the coefficient of the zero order term of the elliptic operator is non-negative, at 
least if a source term is present (see for example [39l Sect. 3]). We give here a 
reformulation, used in several places below, which instead requires non-negativity 
or positivity of the principal eigenvalue. 

Lemma 4.2. Let L be a linear elliptic operator of the form ( f70j) on a compact 
manifold. Let A and <p > be the principal eigenvalue and eigenfunction of L, 
respectively, and let ip be a smooth solution of Lip = f for some smooth function 
f > 0. Then the following holds. 

(i) If X = 0, then f = and ip = C <f) for some constant C . 

(ii) If A > and f ^ 0, then ijj > 0. 

(Hi) If A > and f = 0, then ip = 0. 

Remark. Clearly analogous results hold for / < 0. 

Proof. For a positive principal eigenfunction <p of L, we define \ by ip = x4> and 
an operator T[(f)] by the first equation in 

V[4>) X = -D A (tfh AB D A x) + 2<f?t A D AX + \4> 2 X = <Pf (19) 

while the second equality follows by computation. The strong maximum principle 
[39| Theorem 3.5], applied to ffl9l) gives the desired results. □ 

We end this section with a result which is essentially a straightforward appli- 
cation of linear algebra. 

Lemma 4.3. Let L be a second order elliptic operator on S of the form (f77|). Let 
X,(p,4>^ be the real principal eigenvalue, and the corresponding real eigenf unctions 
of L and its adjoint L'. Let P be the projection operator defined by 




and let Q = I - P. Then 
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(i) L = XF+A, where A has spectrum a (A) such that for some constant cq > 0, 

min | /x — A | > c (20) 

(ii) For any s > 2, p > 1, there is a constant C such that the inequality 

}u\\ W s, P < C\\(A - XQ)u\\ W s-2, P (21) 



holds, where W s,p is the usual Sobolev space of functions with s derivatives 
%n L P (S). 

Remark. We refer to the constant Co as the spectral gap. 

Proof. We shall consider complex eigenvalues and eigenvectors, and use the sesquilin- 
ear L 2 pairing 

(u, v) = I uv. 



JS 

It is straightforward to check that LP = PL = AP, and if we define the operator 
A by A = LQ, so that 

L = AP + A, 

then 

AP = FA = 0, AQ = QA = A, A f P f = 0, A t Q t = A\ (22) 

where we used that P"^ is the projector onto </>' i.e. P'/ = (f>'((f),f)/(4>\(t>). It 
follows that the range of A lies in the space orthogonal to </>'. 

Let ip be an eigenvector for L corresponding to a non-principal eigenvalue 
/i 7^ A. We find that 



^^) = (0 t ,L^) = (Lt0t^) = A (0t ;V ;). 



Since ji ^ A by assumption, we conclude that (0^,^) = 0. Combining the 
previous facts, it follows easily that the spectrum a (A) of A satisfies 

°IA) = ML) \ {A}) U {0}. 

Since A is the principal eigenvalue of L, so that A < Re(/i) for any non-principal 
eigenvalue fi, (J20l) follows from the fact that L has discrete spectrum. 



To prove point (ii) , we first note that by construction, A — AQ is a Fredholm 



operator which maps the space orthogonal to $ into itself. By point (i) , the 
spectrum of A — AQ is bounded from below by cq > 0. Thus, by the Fredholm 
alternative A — AQ : (0^) _L — > (0^)" L is invertible on Sobolev spaces. □ 
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5 Stability of MOTS 



The concept of stability of a marginally outer trapped surface with respect to 
a given slice E is a central issue of our paper and crucial for the application to 
existence of marginally outer trapped tubes. We first briefly comment on stability 
of minimal surfaces embedded in Riemannian manifolds (E, 7), disregarding any 
embeddings in spacetimes. Letting m l be the unit normal to S, the stability 
operator L m (() = b~^ m p where p is the mean curvature of S reads (e.g. [4"T] ) 
L m ( = — [R^rrtm? + KijK l: >)(, where and Rij and are the Ricci tensor 

of (E, 7) and second fundamental form of S respectively. This expression also 
follows from ([5]) by taking S immersed in a time symmetric spacelike hypersurface 
(E, 7) and choosing v = m. As L m is self-adjoint its principal eigenvalue A can 
be represented as the Rayleigh-Ritz formula ( TTTT) . In terms of the latter, the 
standard formula for the variation of the area A of a minimal surface along a 
vector v l = ipm 1 gives 

6 2 V A = 6 2 fr ]s =l' ^8 vPns = I ^L m tfj Vs > \ ! V>V (23) 

The minimal surface is called stable if A > 0, and (123]) together with (TTTT) 
shows that this is equivalent to 5% A > 0, for all smooth variations S v . 

We wish to characterize stable MOTS embedded in arbitrary hypersurfaces in 
spacetime in a similar way as stable minimal surfaces in Riemannian manifolds. In 
the general case we lose the connection between stability and the second variation 
of the area (except if the variation is directed along l a ). Using the discussion 
in the previous section, we can now define stability either in terms of a sign 
condition on suitable first variations of 9 at the MOTS or, in view of (j5J) by 
a sign condition on the principal eigenvalue of the stability operator. In either 
case, we wish to establish a relation between both concepts. Here we choose the 
variational formulation as definition, from which we show the properties of the 
principal eigenvalue. 

We could obtain a variational definition by replacing ( TTTT) by one of ( jl2l) . ( TT3j) 
or ( TT6j) but the resulting definition does not seem very illuminating. Therefore, 
instead of defining stability in terms of sign requirements on integrals overS for all 
positive variations of 9, we now require the existence of at least one variation along 
which 9 has a sign everywhere on S. This definition also enables us to introduce 
an important refinement, namely to distinguish whether this preferred variation 
of 9 is just non-negative everywhere, or even positive somewhere. (We slightly 
expand our earlier presentation P3J. We recall that v a is linearly independent of 
l a everywhere on S. 

Definition 5.1. A marginally outer trapped surface S is called stably outer- 
most with respect to the direction v iff there exists a function ip > 0, ip ^ 0, on 
S such that 5^ v 9 > 0. S is called strictly stably outermost with respect to the 
direction v if, moreover, S^ v 9 7^ somewhere on S. 
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We will omit the phrase "with respect to the direction v" when this is clear 
from the context. We now establish the connection between stability and the 
sign of the principal eigenvalue. 

Proposition 5.1. Let S C S be a MOTS and let X be the principal eigenvalue of 
the corresponding operator L v . Then S is stably outermost iff A > and strictly 
stably outermost iff A > 0. 

Proof. If A > 0(> 0), choose ip in the definition of (strictly) stably outermost as 
a positive eigenfunction corresponding to A. Then 5^9 = L v (p = X(f> > 0(> 0). 
For the converse, we note that from Lemma 14.11 the adjoint L\ (with respect to 
the standard L? inner product ( , ) on S) has the same principal eigenvalue as 
L v , and a positive principal eigenfunction <fy. Thus, for tp as in the definition of 
(strictly) stably outermost, 

A(0t, V>> = (LW, if>) = (0 t , L v ip) > 0(> 0), 

Since (<p',ip) > 0, the lemma follows. □ 

The stability properties discussed above always refer to a MOTS S with re- 
spect to a fixed direction v a normal to S or a fixed hypersurface to which v a is 
tangent. This raises the question as to how the stability properties depend on 
this direction. Using definition ([2]) we can change v a by adjusting the function V, 
which can take any value, and study the dependence of the principal eigenvalue 
A on this function. This yields the following 

Proposition 5.2. Let S be a MOTS and let X v , X v > be the principal eigenvalues 
of the stability operators in the directions v a and v' a defined by with some 
functions V and V respectively. Let <$' and (jr be principal eigenf unctions of L v i 
and Ll respectively. Then 

(X v - A„O(0 t , 0') = (0 f , (V - V)W<ff). (24) 



Proof. From (Q it follows that, 

L v = L v i + {V - V)W (25) 

Hence, = (<ft, (L v — L v > — (V — V)W) 4>') and the proposition follows using 
point (ii) of Lemma 14.11 □ 

Applying some trivial estimates to ff24|) we find 



A„/ + mi[(V' - V)W] <X V < X v > + swp[{V - V)W]. (26) 
5 s 



The inequality fl2"61) implies in particular that, in spacetimes satisfying the null 
energy condition, a MOTS which is stable or strictly stable with respect to v a 
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will have the same property with respect to all directions "tilted away" from the 
null direction l a defining the MOTS, and it puts a limit to the allowed amount of 
tilting of v a "towards" l a which preserves stability or strict stability. In particular, 
if a MOTS S is (strictly) stable with respect to some spacelike direction, then it 
is also (strictly) stable with respect to the null direction — k a complementary to 
l a in the normal basis {l a , k a }. 



6 The graph representation of MOTS 

We now assume that (Ai,g) is foliated by hypersurfaces S t defined as the level 
sets of a smooth function t. Assume also that one of the elements of the foliation, 
say S , contains a MOTS So- We further assume that S is transverse to the 
null normal vector l a on <S, i.e. T p M = T p E © (l a \ P ) for all p E S . However, 
for most results we do not assume any specific causal character for the foliation 
Xi and we allow leaves which are spacelike, timelike, null or which change their 
causal character. The transversality above allows us to fix l a so that l^d^t) = 1. 

The main goal of this paper is to examine under which conditions there is a 
MOTT g such that S t = Q n E t is a MOTS for all t. For this purpose it is useful 
to define an abstract copy S of S detached from spacetime. We are looking for 
a one-parameter family of smooth immersions Xt '■ S — > Tl t for t in some open 
interval I 3 such that S t = Xt(S) is a MOTS and the map 

* : I xS > M 

(t,p) ► i Si o Xt (p), 

is smooth, where i-£ t : S t M. is the natural inclusion. In other words, the col- 
lection of MOTS should depend smoothly on t. Using arbitrary local coordinates 
y % on E t , x A in S, and writing = {t, y 1 }, the immersion ^ takes the local form 

W(t,x A ) = (t,xl(x A )). (27) 
For the null expansion we obtain 

V-Wt K AB -h t l»\l tABdxC dxAQxB K P \ x= * a dxA dxB ) , V*) 

(S) c 

where h t AB is the induced metric on S t , T). AB are the corresponding Christoffel 
symbols and T^ p are the spacetime Christoffel symbols. The vector la(x A ) satisfies 
^a =0 = l a and solves 

^ = 0, 1^ = 0, J'"0„f = l. (29) 

under the condition that l l a depends continuously (and hence smoothly) on t. The 
condition we want to solve is 6 = which is a scalar condition on the immersion 



16 



Xt- Since the problem is diffeomorphism invariant, we need to choose coordinates 
in order to convert the geometric problem into a PDE problem. The idea is to 
construct a suitable coordinate system adapted to the initial MOTS So and to 
restrict the class of allowed MOTS to be suitable graphs in this coordinate sys- 
tem. Since So may have self-intersections, the coordinate system and the graphs 
we will use are only local in the sense that for any point p e S, there is a space- 
time neighbourhood V p of its image such that the spacetime metric takes a 
special form in suitable coordinates adapted to the surface (strictly speaking, to 
the connected component of Sq n V p containing In the following lemma 

we introduce this adapted coordinate system. In order to avoid cumbersome no- 
tation, we will assume S to be embedded for this lemma. Since the applications 
of this result below will always be local on So in the sense just described (and all 
immersions are locally embeddings) this suffices for our purposes. 
The following lemma does not require So to be a MOTS. 

Lemma 6.1. Let (Ai,g), H t and t be as before. There exists a spacetime neigh- 
bourhood V of a smooth, embedded closed spacelike 2-surface So C S C M., local 
coordinates {t, r, x A } on V adapted to the foliation {S t } and functions Z,ip,r] A 
and Uab such that the metric can be written as 

g^dx^dx" = 2e z dtdr + cpdr 2 + h AB (dx A + r] A dr) (dx B + r] B dr) , (30) 

where S H V = {t = 0, r = 0} ; Z(t = 0,r = 0,x A ) = and h^B is a positive 
definite (n — 2) -matrix. The function tp may have any sign (even a changing one) 
reflecting the fact that the hyper surf aces T, t , which coincide with the t = const 
surfaces, are of any causal character. 

Proof. Let x A be a local coordinate system on So and let v a be a smooth, nowhere 
zero vector field on S orthogonal to this surface, tangent to S , and satisfying 
v a l a = 1. We extend v a to a smooth vector field (still denoted by v) in a neigh- 
bourhood of Sq within S and define coordinates (r, x A ) on this neighbourhood 
by solving the ODE (again with slight abuse of notation) 



In some tubular neighbourhood of So within E , this solution defines a smooth 
coordinate system in which we have v a d a = d r . Furthermore, for small enough 
r the sets S ro = {r = r } define spacelike, closed, codimension-two surfaces 
embedded in So and diffeomorphic to <So- After restricting the range of r if 
necessary, we can also extend l a to a nowhere vanishing null vector field on So 
orthogonal to each level surface {S r } and satisfying l a d a t = 1. 

We finally consider the null geodesies with tangent vector l a whose "length" is 
fixed by l^d^t = 1 everywhere in a suitably small neighbourhood of So in M.. On 
this neighbourhood V define functions {r, ar 4 } by solving l a (r) = 0, l a (x A ) = 
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and so that r|s and x A \^ coincide with the functions with the same name defined 
above. The functions {t,r,x A } define a coordinate system on V. Since l a d a = d t 
everywhere we immediately have gu = for the metric. On So = {t = 0} we have 
(d t , 8 x a) = 0, where (, ) denotes scalar product with g. The geodesic equation 
fVal' 3 = UP becomes dtg^t = bg^ t in this coordinate system {b need not be zero 
as the null vector l a has already been chosen). Hence g tx A = and g tr > on 
this neighbourhood (because g tr — 1 on S ). □ 

In terms of the coordinates (130]) . we can consider local graphs St given by 
r = f(x A ) on the surface {r = 0} C E t . Restricting the allowed immersions ff27j) 
to those local graphs, ff28l) becomes an operator on /, which we call 9 t [f]. We 
formulate the ellipticity property of this operator in the subsequent lemma. We 
use the shorthands f A = d x Af and p = \^pf A f A + (l + ^ A /a) ]r=/> where (p and r\ 
are the metric coefficients in (130|) and capital Latin indices are moved with Iiab 
and its inverse h AB '. 

Lemma 6.2. Consider codimension-two immersions defined locally by a smooth 
function f(x A ) according to 



Xt [f]:S^S t cM 



x 



1 _ (t = t ,r = f(x A ),x A = x A ). 



If f fc < C for a suitable positive constant C, the operator 6 ^[f] is quasilinear 
and uniformly elliptic. 

Proof. We give explicitly some steps of this straightforward calculation, which 
requires some care if Ef is not spacelike. We define the quantities 

7 = ^-V z (p-^l + V A f A )-l) , u = p-h z \ r=f , (31) 

\ / r =f 

which are smooth even at points where <p vanishes (i.e. where Ef becomes null). 
The one-form l l a satisfying ( 1291) is 

/* dx a = 1 dt + v (dr - f A dx A ) , 

For the metric induced on the St and its inverse we have 

htAB = (h A B + ffAfB + fAVB + fsVA + ^o/a/b) \ r=f , 

h AB = h AB + p- 1 [f D f DV A V B - V f A f B - (1 + V D f D ) ( v A f B + f A v B )} \ r=} . 

Note that hf B is positive definite wherever p > 0. The operator t [f] can now be 
written explicitly 

9 t [f\ = -vhf B f AB - (ll(T? r h? B )fAfB + K V r A h AB fB + l^ AB h AB )\ r=f , (32) 

where f A B — d x Ad x Bf. Choosing C small enough it follows that p > e for 
a positive e. The assertion of the lemma is verified easily by estimating the 
eigenvalues of vh AB . □ 
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Remark. We note that for quasilinear elliptic equations such as 9 t [f] = there 
hold regularity results. In particular, if we required C 2,a for the function / 
instead of smooth as in the definitions above, the latter differentiability would in 
fact follow (see e.g. Sect. 8.3 of Evans [29J ) . 

We also remark that the stability operator (J5J) coincides with the linearization 
of the quasilinear operator (1321) and can be obtained from the latter by making 
/ infinitesimal. However, the expression (J32l) requires a choice of coordinates. 
On the other hand (jHJ) was derived in a covariant manner and therefore holds 
independently of the coordinates on S or on the spacetime. 

7 Barrier properties of stable MOTS 

Let Q be a marginally outer trapped tube, i.e. a 3-surface foliated by marginally 
outer trapped surfaces. Then there exists a positive variation along Q such that 
<Vu$ = 0. According to the definition in section [5] the MOTS foliating a MOTT 
are stable with respect to the MOTT, but they are not strictly stable as 5^ v 9 = 
L v ip = implies that the principal eigenvalue vanishes. In general situations, 
since the variation of 9 is essentially a derivative of the (9's of adjacent surfaces, we 
expect that embedded strictly stable MOTS are "barriers" , i.e. local boundaries 
separating regions containing weakly outer trapped and weakly outer untrapped 
surfaces. The difficulty in showing this lies again in the fact that strictly stable 
MOTS had to be defined in terms of a sign condition on a single variation, while 
they should be barriers for all weakly outer trapped surfaces in a neighbourhood. 

Below we give a result on this issue, which in fact requires the full machin- 
ery developed in the preceding sections, namely the properties of the principal 
eigenvalue and eigenfunction of the stability operator as well as the maximum 
principle applied to the quasilinear elliptic equation (|32|) representing the MOTS 
as a graph. The maximum principles for non-linear operators normally require 
uniform ellipticity (c.f. eg. [391 Sect. 10]) which is not the case for operators 
of prescribed mean curvature such as ( l32i) when f fc is not small (c.f. Lemma 
16. 2p . For this reason we prove the following theorem "from scratch". To state 
the result, we recall from [1] the definition of "locally outermost". Note that in 
the definition and in the theorem below we require that the MOTS is embedded 
rather than immersed as in the previous sections. 

Definition 7.1. Let S be a MOTS embedded in a hyper surface E. S is called 
locally outermost in E ; iff there exists a two-sided neighbourhood of S in E 
such that its exterior part does not contain any weakly outer trapped surface. 

Theorem 7.1. (%) An embedded, strictly stably outermost surface S is locally 
outermost. Moreover, S has a two-sided neighbourhood U such that no 
weakly outer trapped surface contained in U enters the exterior of S and no 
weakly outer untrapped surface contained in U enters the interior of S. 
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(ii) A locally outermost surface S is stably outermost. 



Proof. The first statement of (i) is in fact contained in the second one. To show 
the latter, let be the positive principal eigenfunction of L v . Since L v (p > 
by assumption, flowing S in S along any extension of <pv a within £ produces a 
family S a , a G (— e, e) for some e > 0. By choosing e small enough, the S a have 
Q\s„ < for o G (— e, 0) and 0\s a > for cr G (0, e). We can now take U to be the 
neighbourhood of S given by U = U ae ^_ e ^)S CT . We first express the expansion of 
the level sets of o in terms of connection coefficients. Setting f = a = const. > 
in (132]) . only the last term survives and we have 

o<e t [o-] = -ilr%h AB \ f=a (33) 

and analogously for a < 0. 

Now let B be a weakly outer trapped surface ( i.e. 6\g < 0) contained in 
14 which enters the exterior part of U. Let p be the point where o~\g achieves a 
maximum value. In a small neighbourhood of p, B is a graph given by a function 
/ which achieves its maximum at p. From ( |32l) . 



O>0\ B = W\ = ~vh? B f AB - (ll(T? r ht B )fAfB + 2l^ rA h AB f B + llV%h AB )\ r=} 

(34) 

At p, /a\ p = and h /ab|p < 0. Since a = const and B have a common 
tangent plane at p, the last term in (JMj) coincides with (J3"3~j) . which yields the 
contradiction (we also use z/L = e z \ 



9 t [f] p = -e z h AB f AB \ p - llT%h AB \ p > 0. (35) 

Hence B cannot enter the exterior part of U, and analogously weakly outer un- 
trapped surfaces cannot enter the interior. 



To show (ii), assume S is locally outermost but not stably outermost. From 
Proposition 15.14 the principal eigenvalue A is then negative. Arguing as above 
one constructs a foliation outside S with leaves which are outer trapped near S, 
contradicting the assumption. □ 

Returning to the example of a MOTT mentioned above, it follows that a 
MOTS within a MOTT is not only not strictly stable with respect to the direction 
tangent to the MOTT, as we already knew, but moreover not locally outermost. 

The barrier arguments suggest that a region of a slice T, t bounded by an 
outer trapped surface Si and outer untrapped surface S2 contains at least one 
MOTS. This has been proven recently by Andersson and Metzger [31], based on 
an argument proposed by Schoen 
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8 Symmetries 



As examples one normally considers spacetimes with symmetries, i.e. invariant 
under the action of some symmetry group. If the symmetry is spacelike, it is 
natural to consider spacelike foliations with the same symmetry. In order to locate 
the MOTS contained in the leaves and to compute the principal eigenf unctions 
and eigenvalues, it is useful to know whether the MOTS and the eigenf unctions 
inherit the symmetries from the ambient geometry. In this section we address 
these questions, starting with the MOTS. 

Theorem 8.1. Let H be a local isometry of (A4,g) (i.e. C^g^ v = for the 
Lie- derivative w.r. to the corresponding Killing field £ a ), and let S be a 
MOTS which is stable with respect to a normal direction v a such that the normal 
component £ ±a of £ a satisfies £ ±a = if)v a for some function if). Then either 5 
leaves the MOTS invariant (i.e. £ a is tangential to the MOTS), or S(«S) is a 
MOTT. 

Proof. Clearly £ a leaves the MOTS invariant iff £ ±a and hence if) are identically 
zero. Assume this is not the case. From ([7]), = 5^9 = 5^ v 9 = L v ip shows that 
if) is an eigenfunction of L v with eigenvalue zero. Since the MOTS is stable with 
respect to v a , if) is the unique (up to a constant) principal eigenfunction of L v and 
has therefore a constant sign, which after reversing £ a if necessary can be taken 
to be positive. This implies that S in fact generates a MOTT as stated. □ 

Note that this proof shows in particular that if S is strictly stable then S 
must remain invariant under the isometry. 

If S lies in some hypersurface S invariant under the isometry, then Theorem 
18.11 implies that either H leaves S invariant, or H(«S) C S. Clearly, in the second 
case S is not strictly stable. Moreover, it is not locally outermost in the sense of 
Definition 17.11 

We also remark that the range of validity of Theorem 18.11 can be extended 
with the help of Proposition 15.21 It suffices to require that the MOTS is stable 
with respect to any direction normal to S which lies "between Z M and £ ±fl " (which 
is in fact a conical segment). 

The following theorem on the symmetry of the principal eigenfunction ( "ground 
state") is well-known for self-adjoint operators, with numerous physical applica- 
tions. We have formulated it here for the general linear elliptic operator (fTU|) . 
and therefore it holds in particular for the stability operator (J5]). 

Theorem 8.2. If L is invariant under a 1-parameter group of isometries gener- 
ated by r] a , the principal eigenfunction <f> is invariant as well. 

Proof. The invariance of L under the isometry implies that this operator com- 
mutes with the action of the corresponding Lie derivative C v . Hence, if <f> is an 



21 



eigenfunction with principal eigenvalue A, it follows that 

LC v (f) = C v L(f) = XC v (j) (36) 

which means that C n <f) is a principal eigenfunction as well. But from Lemma 
I4.1( i). the latter is unique up to a factor, i.e. C^cj) = a<p for some real number 
a. This expression integrates to zero on S because C v (p = Da(<Pti a ). Since has 
constant sign, it follows that a = 0. This proves the assertion. □ 

As an example we consider axially symmetric data on E with Killing vector 
n a = d/dip, which has in addition a (t, tp) symmetry, i.e. invariance under simul- 
taneous sign reversal of tp and t. This is the case, in particular, for data on a 
t = const, slice of the Kerr metric in Boyer-Lindquist coordinates. It turns out 
that the vector sa introduced in Sect. [3] is proportional to the axial Killing vector 
and divergence- free, i.e. sa = za in the notation of Sect. HI By Theorem 18.21 
the principal eigenfunction is axially symmetric and therefore satisfies a second 
order ODE. Finding the principal eigenvalue and eigenfunction is therefore equiv- 
alent to solving a one-dimensional Sturm-Liouville problem. Moreover, the first 
order term in the stability operator vanishes when acting on this eigenfunction. 
However, this does not imply that the stability operator is self-adjoint in this 
situation or that the principal eigenvalue coincides with any of the symmetrized 
eigenvalues \ s or \ z . 

9 Existence and properties of MOTTs 

The main objective of this paper is to show that MOTS "propagate" from a given 
slice to adjacent leaves of a given foliation to form a MOTT. Setting 6 = in (1251) 
determines the location of a MOTT in terms of immersions. These may be viewed 
as defining a family of quasilinear elliptic equations which do not contain any 
derivatives along the presumptive "evolution" direction. Since we have assumed 
that iS is marginally outer trapped, the elliptic equation is satisfied for t = 
and we can adopt a perturbational approach. In section O we have derived the 
"graph" representation (1321) in a special coordinate system, which will be used 
in the existence proof. The proof makes use of a general existence result [32] 
which shows the existence of solutions near a known solution of a one-parameter- 
family of non-linear differential equations (or systems) of arbitrary order and 
arbitrary type whose linearization is elliptic. Since we will apply this result to 
a 1-parameter family of quasilinear elliptic equations for which regularity results 
hold, we will formulate the result here for smooth functions only. Equations of 
general type require some care regarding differentiability. 

Lemma 9.1 ([321 [33]). Let 

F{x,u,du,d 2 u;t) = (37) 
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be a 1-parameter family of quasilinear, second- order differential equations, where 
F is a smooth function of all arguments. Assume that Uq is a smooth solution of 



(31) for t = to, and that the linearized equation around uq, Lv = f is elliptic and 



has a unique solution for any smooth f . Then (31) has a unique smooth solution 
for t E (to — e, to + e) for some e > 0. 

This result is an easy application of the implicit function theorem; alterna- 
tively, the latter theorem can also be applied directly to (I32l) . as we suggested in 
PQ. The result is, in view of the above discussion, that a smooth horizon exists in 
some open neighbourhood of the given slice. More precisely, we have the following 
theorem and a corollary. 

Theorem 9.2. Let (M.,g a p) be a spacetime foliated by smooth hypersurfaces Y, t , 
t G [0, T], and assume that So contains a smooth, immersed, strictly stable MOTS 
So- Then for some r G (0,T] there is a smooth adapted MOTT 

g [0 , T) = $(<S x [0,r)) 

such that for each t G [0,r) ; S t = Q(So,t) is a smooth, immersed, strictly stable 
MOTS with 5j C E t . 

Corollary 9.3. The MOTT Q[o, T ) constructed in Theorem \9.S\ is nowhere tangent 
to the leaves S t of the foliation. 

Proof. Let {t, r, x A } be the local coordinate system introduced in Lemma [67T1 and 
restrict the allowed MOTT to be local graphs on S as described in section [6j 
Thus, we are looking for functions f(t,x A ) which satisfy 0t{f] = for all t near 
t = and satisfying f(0,x A ) = 0. In order to apply Lemma 19.11 we only need 
to check that the linearization of 9 t [f] is elliptic and invertible. The immersion 
$ defining the MOTT is given in this coordinate system by (t,r = f(t,x A ),x A ). 
Recall that l a d a = d t and that on So the vector v a d a = d r satisfies l a v a = 1. 
Linearizing the operator around (t = 0, / = 0) corresponds to fixing t = and 
making / infinitesimal, or more precisely to evaluating L(f') = d e 9 t =o[fe]\e=o 
where f t depends smoothly on e and satisfies f t \e=o — 0, d e f t \ e=0 = f. It follows 
that L corresponds to performing a geometric variation of 9 along the vector f'd r . 
The general variation formula in section [3] gives L(f') = 6f v = L v (f). Since by 
assumption, So is strictly stably outermost with respect to v, it follows that L v 
and therefore L, is invertible. Existence of f(t,x A ) for small t solving 9t[f] = 
follows readily from Lemma I9TT1 Moreover, the "evolution" vector to the MOTT, 
q a = $*(<9f) = dt + f'd r is by construction nowhere tangent to £f, which proves 
the corollary. □ 

As a complement to this Corollary we now discuss the "tangency" property 
of the MOTT when the principal eigenvalue A T goes to zero, assuming smooth 
convergence of the S t to a stable MOTS S T . 
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Theorem 9.4. Let (M,g a /3) be a spacetime containing a smooth reference fo- 
liation {Ejjjg^T]- Assume that S contains a smooth, immersed, strictly stable 
MOTS Sq. Assume furthermore that the adapted MOTT G[o, T ) through So con- 
structed in Theorem \9.B is such that as t — ► r, the surfaces S t converge to a 
smooth, compact, stable MOTS S T . Let X t be the principal eigenvalue of the sta- 
bility operator of S t and <p\ the principal eigenfunction of its adjoint. Assume 
that A i ~ 1 (0j, W|s t ) has a limit (finite or infinite) as t — > r. 

Then the closure G[o, T ] — ^[o,r)U5 T is an adapted MOTT. If in addition X T = 
and ((fil, W\s T ) 7^ 0, then G[o,r] is tangent to S r everywhere on S T . 

Corollary 9.5. If the null energy condition holds on S T , X T = and W\$ T ^ 
somewhere, then Q[o lT ] is tangent to S r everywhere on S r . 

Proof. If A T > 0, then by uniqueness the MOTT through S T constructed using 
Theorem 19.21 agrees with G[o, T ) f° r t G [0, r), and hence the closure G[o, T ] is an 
adapted MOTT. It remains to consider the case where A T = 0. 

By construction, S t — E t fl G[o, T ) is a MOTS for t G [0, r). Let v a be the 
normal to St tangent to £t (we drop the subindex t for simplicity). Consider 
the "evolution" vector q a of the MOTS within G[o,t), given by q a d a = ty*(dt) = 
dt + (dtf)d r . Recalling that l a d a = dt it follows that the normal component of 
q a can be decomposed as 

q ± a = r + uv a ^ 

for some function u. The variation of 9 along q a must vanish (this is precisely 
the condition that q a is tangent to the MOTT). Hence, u must satisfy 

6 g 6 = 5 l+uv 9 = -W + L v u = 0, (39) 

c.f. Lemma [3. 2[ 

Consider the equation L v u = W. Referring to Lemma H~3l we decompose u = 
Fu + Qu and write this as u = zip + w for some real z. We have L v u = Xzcf) + Aw, 
where A = L V Q as in Lemma 14.31 It follows that 

—IS 

and 

Aw = QW. 

Since for t G [0,r], L v is uniformly elliptic and its coefficients are uniformly 
smooth, [121 Theorem 3.1, p. 208] applies to show that the spectral gap is lower 
semi-continuous in t, i.e. liminf co(t) > cq(t). At r, the spectral gap is positive. 
Thus, it follows that there exists e > and a constant Co > such that the 
spectral gap estimate (|2"U1) holds for t G (r — e, r}. Recall that by construction 
Qw = w. Since by assumption A > for t < r, we may apply (|2~T|) with A = 0, 
to conclude that Qw, and hence also w, is uniformly bounded in Sobolev spaces. 
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Suppose now that (0t, W\ Sr ) ^ and A T = 0. Then it follows that z — > oo, 
as t — > r, and hence that u diverges. We see further that in this case, 

limz _1 M = (p T > 0, 

t^T 

where <p T is the eigenfunction of A T . Define a normalized evolution vector q a = 
z~ x q a . Then 

q a d a = z~ x d t + z~ l ud r . 

It follows from the above that as t /* r, the vectors g a converge smoothly to 
a limit, which is proportional to d r . By construction q a is a smooth vectorfield 
tangent to the MOTT Q[o )T ), and we have now shown that q a extends smoothly 
to the closure G\o, T ] = G[o, T ) U St- In order to demonstrate that the closure is 
indeed a MOTT, it remains to renormalize the time parameter. 

Thus, let 4> s : <So — > £/[o,r] be defined by the flow of the vector field q a , and let 
/ = [0, s*] denote the parameter interval required for this flow to sweep out G[o, T \- 
This construction defines a smooth monotone map a : I — > [0,r]. We can now 
extend the immersion defining the MOTT £/[o )T ) to an immersion $ : Sq x / — > M. 
such that for s G /, we have «S CT ( S ) = $(«Sq, s) C Scr( s ) for s £ I. It is clear from 
the construction that this defines a MOTT £/[o,t] = @[o,r) U »5t- 

It remains to consider the case when W|5 T ) = 0. In this case, by applying 
the same argument as above, we find that u may, or may not diverge as t — ► r, 
depending on the detailed behavior of A and W as t — ► r. In particular, we see 
from (jlDj) that if A - ^^''", W) diverges as t — >■ r, then also 2; diverges, and we are 
essentially in the situation considered above. On the other hand, if X~ l {(f)\W) 
tends to a bounded limit as t — > r, then we are in a situation which is analogous to 
the case with A r > 0. In either case, the normalized evolution vector q a converges 
as t — > t, and an argument along the lines above shows that G[o, T ] — G[o,t) U $t is 
an adapted MOTT. □ 

In [I] we claimed that the evolution of the MOTS can be continued "as long 
as the MOTS remain strictly stably outermost". We emphasize, however, that 
in p] the MOTS were taken to be smooth and embedded by definition. Moreover, 
we have tacitly assumed that the MOTT does not "run off to infinity" . Hence, 
our control of the evolution may end not only when A goes to zero but also when 
the MOTS develop self-intersections or when we lose compactness or smoothness. 
In the present setup, we do not worry about MOTS developing self-intersections 
as we have allowed them from the outset. However, to show in particular that 
existence continues up to and including E T (which was assumed in Theorems 19.21 
and 19.41) we have to exclude the other pathologies. 

To avoid that the MOTT "runs off to infinity" we simply require that it is 
contained in a compact subset of Ai. Note that it may still happen that the 
area of the MOTS St grows without bound as we approach «S T . In view of the 
curvature estimates of [34J, this can happen in the four-dimensional case only if 



25 



the MOTS "folds up" sufficiently. Assuming a uniform area bound, we have the 
following result, which follows from the work in [34] . 

Proposition 9.6. Let (Ai,g a /3) be a spacetime of dimension 4, with a spacelike 
reference foliation {Et} fe [ 0jT ] . Assume that E contains a smooth, immersed, 
strictly stable MOTS So. Let Q[o, T ) be the adapted MOTT through So constructed 
in Theorem \9.Sl and let St = G[o, T ) H St be the leaf of G[o, T ) St. Assume that 
for t G [0, r), the leaves S t have uniformly bounded area, and are contained in 
a compact subset of M. Then the S t converge as t — > r to a smooth compact 
surface S T which is a smooth, immersed, stable MOTS. 

If the dominant energy condition is satisfied, the area of the limit set of the 
MOTT in fact stays bounded as long as the MOTS stay strictly stable, as a 
consequence of the following lemma. (This computation is known in the context 
of the topological results [2D])- Recall that v a is everywhere linearly independent 
of P. 

Lemma 9.7. In a 4-dimensional spacetime (M.,g a p) i> n which the dominant en- 
ergy condition holds, let S be a MOTS which is strictly stable with respect to a 
spacelike or null direction v, with principal eigenvalue A and area \S\. Then S is 
topologically a sphere and \\S\ < 47r. Moreover, if \\S\ = Att then S has constant 
curvature, i.e. Rab = A/iab- 

Proof. We take to be the principal eigenfunction of L v , and we call y A = 
—<P~ 1 Da(P + s A and y 2 = yAV A - Taking if) as the eigenfunction in fl^D, we obtain 

A = D A y A - \y\ 2 + l -R s - Y, (41) 

where Y has been defined in fl3j). Integrating (14ip and using the Gauss-Bonnet 
theorem gives 

X\S\ = 2n X - Jy + Y), (42) 

where \ is the Euler number. The first assertion of the lemma now follows since 
A > and Y > 0. If A|<S| = 4vr, flU implies y A = and Y = 0. Putting this 
back into (f4T|) . we have Rs = 2 A which implies the statement of the lemma as 
Rab = \RhAB in 2 dimensions. □ 

We remark that the same method shows that in the case A = 0, the MOTS 
can be a torus, which necessarily must be flat (Rab = 0). 

With the help of Lemma [9T71 we can now sharpen Proposition [91)] to formulate 
an existence result as follows. 

Theorem 9.8. Let (M.,g a p) be a spacetime of dimension 4, in which the dom- 
inant energy condition holds and which is foliated by smooth spacelike hypersur- 
faces St, t G [0, T]. Assume that S contains a smooth, immersed, strictly stable 
MOTS Sq. Assume further that the MOTT Q\o T ) through So constructed in The- 
orem \9.^ is contained in a compact subset of Ai and that either 



26 



(i) liminfi^ r X t > 0, or 

(ii) lim^ T At = and limsup 4 ^ T |«St| < oo. 

Then, there is a smooth, compact, strictly stable MOTS S T in S T such that Q[o, T ] — 
G[o, T ) U <$t is an adapted MOTT, which is the closure ofG[o )T )- 

We remark that, as the surface S T C S T itself satisfies the requirements of the 
theorem, the evolution in fact continues in [0, T] "as long as the MOTS St C S f 
stay strictly stable" . 

We now reproduce and slightly extend the result [TJ on the causal structure 
of MOTT foliated by strictly stably outermost MOTS, which must be either null 
everywhere or spacelike everywhere. 

Theorem 9.9. Let (A4,g a p) be a spacetime in which the null energy condition 
holds, which is foliated by smooth, locally achronal (i.e. spacelike or null at each 
point) hypersurfaces E t nowhere tangent to l** , and that S contains a strictly 
stable MOTS So. Then, the following holds for any MOTT Q through S§. 

(i) The MOTT Q is achronal in a neighbourhood of So. 

(ii) IfW does not vanish identically, Q is spacelike everywhere near So- 
(Hi) IfW vanishes identically on So, Q is null everywhere on So- 
Proof. Recall equation ( |39l) for the normal variation vector q ±a = l a + uv a . 



Applying points (ii) and (iii) of the maximum principle, Lemma 14.21 proves the 



results. □ 

We finally comment on a possible alternative construction for MOTTs. Con- 
sider a weakly outer trapped surface S on some initial slice and take the null cone 
emanating from it. By the Raychaudhuri equation ([6]) if W > 0, the null cone 
cuts each subsequent slice on an outer trapped surface which gives a trapped 
"barrier" «Si. If we also assume an outer untrapped "barrier" S 2 outside <Si, 
which in particular always exists near spacelike or null infinity in asymptotically 
flat spacetimes, there is a MOTS in the region bounded by Si and S 2 by the 
results of [101 [31], c.f. the discussion in section [71 In this way, one can show 
the existence of an outermost MOTS on subsequent slices. In general however, 
this collection of MOTS need not be a MOTT because it may jump, as it always 
tracks the outermost MOTS on each hypersurface S t . 
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A Proof of Lemma 3.1 



Due to the linearity of the variation and the fact that g" generates a diffeomor- 
phism of S, which implies 5 q \\9 = q^(9), we can assume q = q ± without loss of 
generality. Since the variation is a local calculation we may assume that all S a 
are embedded. Let k a be a null, future directed normal vector to S a satisfying 
k^l^g a p = —2. We wish to extend Z° and k%, for each value of a to a one pa- 
rameter family of vector fields defined on a neighbourhood of S. Let Q a be the 
null hypersurface generated by k" by affinely parametrized geodesies. Extend 
first 1° to Q a by parallel transport along k%. Then extend 1° away from Q a by 
affinely parametrized null geodesies and finally extend k% by parallel transport 
along I". Notice that, in general, the two planes orthogonal to {I", k"} at each 
point do not define two-surfaces. On S a however they obviously do. Being I" 
a geodesic field, the null expansion 9 a can be rewritten as 6 a {p) = (V a Z")|$(p i(T ) 
for any p G S, where a) is the variation of Sect. [3j Defining U a = d a l™\ a= o 
(with partial derivative taken at a fixed spacetime point), we obtain directly from 
the definition of the variation 

M = ^ a U a + q P VpV a r\ s = V a U a + q (3 V a Vpl a -bG aP l a l f3 + -R aP k a l l3 \ s (43) 

where l a = l® =0 and k a = k® =0 and the Ricci identity has been used in the 
second equality. Let us next determine the divergence of U a . From that fact 
that is null for all a, it follows U a l a = and hence U a = al a + XJ A e a A for a 
as in the statement of the lemma and for suitable functions U A . Here e a A is a 
basis of the orthogonal subspace to I and k at each point. Using the projector 
h a /3 = g a /3 + \{lak(3 + k a lp) it is easily found that the divergence of any vector of 
the form F A e a A is 

V Q (F A e a A ) |s = D A F% (44) 

where we used the fact that l a V ' a k^ = and k a V J^ls = 0. We need to determine 
U A on S: in local coordinates y a (x A , a) for the map $, orthogonality of 1% to S a 
means 

Its cr-partial derivative at a = gives, (V eA q,l) + Ua\s — 0, i.e. Ua\s — 
—D A u + us a- From £V«# = it follows U a V a l p + l a V a IJP = 0, which after 
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multiplication with hp and the fact k a is parallel along l a gives l a d a a = —U A s A - 
Thus, using 



V a U a \ s = -A s u + 2s A D A u + u (D A s A - s A s A ) +a6\ s . (45) 

We next consider the second term in (1431) . i.e. bl^V Q V pl a — |A^V a V pl a '. An 
integration by parts and using that l a is geodesic implies 

l p V a V p l a \ s = -VJ p V p l a \ s = -K%K uA %l v \ s . 

Decomposing k^V pl a = Ql a + r A e a A and using the fact that Q\g = r A \$ = 0, 
another integration by parts gives k^V Q V /3l a \s — l^d^Q — K AB K uAB l^k u \s- In 
order to determine l a d a Q, we first note that — 2Q = /c Q /c /3 V a //3. Taking co- 
variant derivative along l a and using the Ricci identity we find 2Z M <9 M <5|.s — 
l a kH^k v R a ^ v \s- Collecting terms we get 

q a V a V^ a \s = -bK^ B K" AB l,l u + | \K>X B K l ' AB l lx k lJ - h a Wk»R a ^ 

'?46) 

For the last term in (|43p . the definition of the projector h af3 gives R + 2R a pl a k 13 = 
h a ^h^R a ^ v +\l a m p k v R a ^ y . Making use of the Gauss identity h af3 h^R a ^ u \ s = 
R s - H 2 + K^abK^ we get 

R a pl a k \ s = -G a pl a k? + R S -H 2 - K^JC^K + \l a ^fk v R a ^ v . 

1 s 

(47) 

Inserting (]4"5l) . (H6l) and ( 1471) into ( 1431) completes the proof. □ 



B Proof of Lemma 4.1 (sketch) 

The Krein-Rutman theorem states that on a Banach space B, a compact linear 
operator E that maps any non-zero element of a closed cone K (i.e. a topologi- 
cally closed subset of B closed under addition and multiplication by non-negative 
scalars) into its topological interior, necessarily has a unique eigenvector u in the 
interior of K of unit norm. Moreover, the corresponding eigenvalue a is real 
and positive and any other element (3 of the spectrum of E (complex in general) 
satisfies a > |/3| where | • | denotes the complex norm. A proof of this theorem 
can be found in |35j . 

In the case of the elliptic operator L fflOl) . let 5 be a constant satisfying 
5 > sup 5 c and define the operator L' = L + 5. The zero order term is therefore 
positive everywhere and the PDE L'f = g admits a unique solution in C 2,a (S) 
for any g G C°' a (S) = B. Let Q : B -> B be defined by Q(g) = f and let K 
be the set of non- negative functions (which is obviously a cone). The maximum 
principle implies that if g G K and non-identically zero, then / = Q(g) is strictly 



29 



positive everywhere. Thus, all the conditions of the Krein-Rutman theorem are 
fulfilled and there exists a unique non- negative function of unit C 0,a (S) norm 
satisfying Q(<f>) = a<p. Since a is positive and (f) is in the image of Q it follows 
that <p is strictly positive and in C 2 ' a (S). Elliptic regularity implies that <ft is m 
fact smooth. It follows that L<p = (a -1 — 5)<ft, so we have a positive eigenfunction 
(unique up to rescaling) and a real eigenvalue A = a" 1 — 5. 

It only remains to show that any other eigenvalue of L has larger or equal real 
part. This does not follow directly from the Krein-Rutman theorem. However, 
we use the following argument, which we adapt from Evans [29]. Let ^ be a 
(possibly complex) eigenfunction of L with eigenvalue /i. Define u = A 
direct computation gives 

- D A D A u + 2t' A D A u + (A - /i)u = 0, (48) 

where t' A = t A — D A (p. Using also the complex conjugate of (l4"8l) a short calcu- 
lation gives 

(-D A D A + 2t' A D A ) \u\ 2 = 2 (Re(/i) - A) \u\ 2 - D A uD A u < 2 (Re(/i) - A) \u\ 2 

Thus, if Re(/i) < A the right-hand side is non-positive and the maximum principle 
would imply \u\ 2 = 0. Thus, Re(/z) > A as claimed. 

Finally, the result on the adjoint is a trivial consequence of the positivity of 
the principal eigenfunctions, c.f. [25]. Explicitly, if and A^ are the principal 
eigenfunction and eigenvalue of L\ it follows 

O=(LVt ; 0)_(0t )jL0) = ( A t_ A ) (0 t ;0) 
and A^ = A as positive functions cannot be L 2 orthogonal. 
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